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This report examines the nature of coherent reflections of radio waves, at near-grazing 
incidence, from a horizontal atmospheric layer on which is superimposed a slight wave 
motion. The existence of reflections of this type is merely postulated, and the development 
then proceeds to examine the consequences of such a postulate with reference to measure- 
ments obtainable in transhorizon propagation experiments. The properties of angle of arrival, 
signal level, fading rate, and Doppler shift are examined, together with their rates of change 
with time. 



1. Introduction 

Certain rapid beam-swinging experiments in trans- 
horizon microwave propagation have been reported 
[Waterman, 1958] as showing an apparently system- 
atic motion of the scattering or reflecting source 
with velocities in excess of wind speed. The sug- 
gestion that these motions might be associated with 
waves 2 on a surface of discontinuity in refractive 
index [Waterman, 1959] has been examined by 
Gossard [1961; 1962] who concluded, on meteoro- 
logical grounds, that the wave velocities were 
incompatible with the angles, slopes, and refractive- 
index discontinuities required. 

This report examines the nature of such hypothe- 
sized reflections with regard to the angle of arrival, 
signal strength, Doppler shift, and fading rates of 
the received signal. It is shown that the location 
of a reflecting facet on a wavy layer does not move 
with the phase velocity of the wave; in fact, it can 
move in the opposite direction, and can have infinite 
velocity. From a single sinusoidal layer and at a 
given instant of time, there may exist one, three, or 
more distinct rays leading to multipath phenomena. 
The variation in signal strength, for the case of a 
single ray, is caused by changes in the angle of 
incidence and by focusing of the energy by the 
layer 's curvature. The Doppler is caused by the 
motion of the reflecting point, both in elevation and 
in azimuth. The effects on the signal strength of 
multiple layers at different elevations are considered. 
For a number of these structures, various signal 
statistics — distribution, power spectrum, time auto- 
correlation — are computed. The main purpose is 
to evolve a theoretical basis for comparing predicted 
results of radio reflections from layers with experi- 
mental measurements. 



» Prepared under Signal Corps Contract DA 36-039 SC-87300. 
2 With regard to the existence of such waves, see- [Gossard 1954, Smyth, 1947 
and du Castel, 1961]. 



2. Geometry of Layer Reflections 

In most transhorizon propagation paths the scat- 
tering or reflecting regions occur near the center of 
the path. If off-path reflections from a layer occur, 
the layer must have an appropriately oriented slope. 
For simplicity of analysis, we consider waves moving 
at right angles to the path and reflection from a 
facet on the midplane between transmitter and 
receiver. The facet's position is then specified by 
an azimuth a, and an elevation e, as seen by the 
receiver. In order that a ray from the transmitter 
be specularly reflected in the direction of the receiver, 
the slope of the facet must be given by 



tan S r = 



tan a • cos e 

shl ( e+ m) 



e + 
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(1) 



in which D is the path length and R is the earth's 
radius. The approximation is valid since the angles 
are small. This is the slope necessary for reflection, 
so designated by the subscript r. It must be equated 
to an expression for the slope of the layer. Consider 
a sinusoidal wave with vertical displacement z, 
amplitude A, and wavelength L, traveling in the x 
direction, normal to the path. It may be represented 

by 

3 =H+A cos (m-^t) (2) 

and its slope is equal to 

, „ dz F2zA . ( 2m; XT 



, 2tA 

H — F~ Sln 



(«-?•) 



(3) 
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if the lateral distance x is replaced by Da/2. (The 
subscript I indicates that this is the layer slope which 
must be equated to tan S r above.) Therefore, at 
the reflecting point, we must have 



2ttA 



e+ 2B 



=—t— sin 



-O-tt) 



(4) 



For a given path, specified by D, and a given wavy 
layer, specified by A, L, £2, and e, a solution of this 
equation for a will determine the azimuth from 
which reflections may be obtained. As time t varies, 
the waves move along the layer and the value of 
azimuth a varies. 

The nature of this variation may be seen by 
solving (4) graphically. Figure 1 illustrates the 
situation in two sample cases. The slope required 
for reflection varies linearly with azimuth — eq (1) — 
and is indicated by the oblique straight line in the 
upper diagram. The slope of a relatively flat sinu- 
soidal layer is indicated by the dashed line. Where 
these two lines intersect, the azimuth of the re- 
flecting point is specified — i.e., (4) is satisfied. As 
this flat sinusoidal layer moves across the path, the 
position of the intersection varies, but there is never 

SLOPE REQUIRED FOR REFLECTION 
-SLOPE OF WAVE 





Figure 1. Graphical method of solution for azimuthal position 
of reflecting point from a sinusoidal wavy layer. 



more than one intersection and consequently never 
more than one reflection. The azimuth of the re- 
flection wanders back and fortli in a distorted 
sinusoidal motion. 

A more pronounced sinusoid, having steeper slopes, 
is indicated by the solid wave in the upper diagram 
of figure 1. Three positions are shown as it moves 
across the path. In one of these, there are three 
intersections with the oblique straight line and con- 
sequently three reflecting points on one layer. The 
lower part of the figure has the same azimuth scale, 
but time runs vertically downward. The dashed 
and solid curves indicate the motion of the reflecting 
point as seen by the receiver for the two cases of a 
flat and steep sinusodial wave. For the latter wave, 
the time interval during which three reflections are 
possible is evident: as one reflection moves off to 
the side in azimuth another appears and immediately 
splits into two reflections, one of which swings 
across the great-circle bearing (zero azimuth) 
and joins the first, where they both disappear, 
while the other swings out to a maximum azimuth 
and then returns more slowly. 

In this figure the rate of change of azimuth with 
time is given by the slope with respect to the vertical 
time scale. Where the curves have a vertical tan- 
gent, the azimuthal velocity is zero; where hori- 
zontal, infinite. Thus the coming and going of the 
multiple reflections are accompanied by momen- 
tarily infinite velocities. 

It is convenient to classify wavy layers in accord- 
ance with the magnitude of the maximum facet 
velocity. For the cases considered here (small 
sinusoidal waves on horizontal layers and moving 
across the path), a parameter rj may be defined 
whose properties are: (1) it is zero for a completely 
flat layer, (2) it equals 0.5 for a layer whose maximum 
facet velocity equals the wave velocity, and (3) it 
equals unity for a layer of just sufficient slope to 
support infinite facet velocities. This parameter is 
related to the layer and path geometry by 



2tA 
= L 



tD ( D\ 
T ' V+2R) 



(5) 



We 



= 0: 
=0.5: 



wave 



facet velocity less than wave velocity 
maximum facet velocity equals 
velocity 
77= 1.0: maximum facet velocity infinite; three or 
more reflections per layer, 
now define waves on a layer as small if 
7/<^0.5, and as large if 77 > 1.0. 

An indication of the variations in the angular 
velocity of the reflecting facet can be seen from 
figure 2, in which azimuth is plotted as a function of 
time for three pairs of different geometrical cases. In 
these cases the amplitude and wavelength of each 
layer have been held constant, while the distance has 
been doubled between Case I and Case II, and again 
between Case II and Case III. In each case, two 
layer elevations are considered, the lower being 
represented by the dashed curve. 

Examining Case I we see that at both elevations 
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Figure 2. Angle of arrival for rays reflected from a wavy layer of form II + A cos ( 12 1 —-j— )• 



the facet velocity never becomes infinite. For these 
two curves the parameter rj is approximately 0.4 and 
0.8 for the elevation angles of 0.5 (leg and 1.5 deg 
respectively. This means that for the lower eleva- 
tion the maximum facet velocity is approximately 
equal to the wave velocity. However, it is in the 
opposite direction; the wave is moving from negative 
to positive azimuths, while the maxim inn rate of 
change of facet azimuth occurs in going from positive 
to negative (at £=50 sec, in the figure). For the 
higher elevation we have a much steeper slope in the 
region near zero azimuth corresponding to a high 
angular velocity of the received ray. 

In Cases II and III, with the increased distances 
rj has become larger than 1.0. As a result we have 
some momentarily infinite velocities and more than 
one reflecting- point. The second case shows as 
many as three received rays during part of the cycle, 
and the third case shows as many as five with never 
fewer than three. 

Figure 2 illustrates the effects that changing the 
distance between transmitter and receiver can have 
on the characteristics of the received signal from a 
given layer. As the distance is increased, the 
angular velocity of the reflecting facet will increase 
until its peak value equals that of the wave; as 
distance is increased further, the peak facet velocity 
w411 eventually become infinite and more than one 
ray will be received. Similar effects would be 
observed for a progressive increase in layer elevation 
or in wave amplitude, or for a decrease in wave- 



length. 



3. Signal Strength 



While there are several factors affecting signal 
strength, we discuss here only some elementary con- 
siderations directly related to the hypothesis of layer 
reflections. First is the magnitude of the reflection 
coefficient for a plane (nonwavy) layer. If the 
"layer" consists of an abrupt change in refractive 
index An from a uniform value below to a different 



uniform value above, the power reflection coefficient 



is 



<&y 



(6) 



where 6 is twice the angle between the plane of the 
interface and the wave normal — i.e., the total devia- 
tion angle of the reflected wave. Compulations 
have been made of the reflection coefficient [Bauer, 
1956] for cases in which the discontinuity is not 
abrupt but is distributed over a distance that is not 
a small fraction of a wavelength. In many reason- 
able cases these computed reflection coefficients do 
not differ markedly from that in (0), and so it will 
be used for purposes of discussion. The magnitude 
of the signal reflected from a layer, then, will vary 
with layer height, other quantities remaining unal- 
tered. Indeed it is interesting to note that the 
power reflection coefficient varies inversely as the 
fourth power of the deviation angle 6. This de- 
pendence is similar to that predicted from the quite 
different model of turbulent blobs [Booker and 
Gordon, 1950]. As regards magnitude, a disconti- 
nuity of one N unit in refractive index (one part in 
10 6 ) and a deviation angle of 2 deg give rise to a 
reflected wave about 55 db below the incident, a 
figure which is comparable with both turbulent- 
model predictions and experimental observations. 

Another effect on signal strength must also be 
taken into account: the focusing factor. Since the 
layer is curved, the reflected rays may diverge less 
strongly or more strongly than in the absence of 
layer curvature. For instance, consider the layer 
previously mentioned. At time t=0, reflection takes 
place on the great-circle plane, and since the layer 
is concave downward, the reflected wave is focused 
toward the receiver, resulting in a higher signal 
strength than for a flat layer with the same discon- 
tinuity. Defining the focusing factor as a ratio of 
the received power after reflection from a curved 
layer to the received power after reflection from a 
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Figure 3. Angle of arrival and signal strength as functions of time for sinusoidal wavy layers with a single reflecting point. 



flat layer, we have 



F= 



1 



l+a-C 



(?) 



where a is the grazing angle and C is the curvature 
of the layer. 

The combined effect of the reflection coefficient 
and the focusing factor leads to signal strength 
curves as seen in figure 3, in which field strength, on 
an arbitrary linear scale, is plotted against time. 
The first row of curves gives the angle of arrival of 
the received ray as a function of time for three 
different cases. The two layers used as examples 
in columns a and b are the same two as in Case I in 
the previous figure; they differ in elevation only. 
The signal strength for the second case is in general 
weaker since the elevation angle, and hence the 
reflection angle, has been increased. During those 
moments of time in which the received ray is rapidly 
crossing the great-circle bearing, there is a peaking 
in the received signal, resulting from reflection of the 
radio wave at the crest of the sinusoidal layer where 
it has maximum downward curvature. During other 
times, the effect of the layer's curvature is less pro- 
nounced. The example in column c shows an alter- 
native comparison with that in column a. In this 
case, the effective layer waviness, indicated by the 
t] parameter, is increased by a near doubling of the 
wave amplitude A, rather than by an increase in 
elevation — all other quantities being held constant. 
The average signal strength is about the same as in 



column a, but the peaking associated with passage 
of the wave crest across the great-circle path is 
greatly intensified, amounting to some 12 db. 

4. Doppler Shift 

The motion of the reflecting facet results in a 
Doppler shift arising from the change in path length 
between the receiver and transmitter. This shift 
in frequency, for a facet on the path-bisecting plane, 
is related to the path geometry and angular velocities 

by 



A 4 D [7 _l_ D \ ^ e _l_ & °1 



(8) 



Once the azimuth versus time curve has been found 
for any particular path geometry, the Doppler may 
be obtained by graphical or numerical techniques. 
Figure 4 illustrates the Doppler for the same three 
cases used in figure 3. The upper row is a repeat, 
to facilitate comparison, of the upper row in figure 3 ; 
columns b and c show the effect of increasing eleva- 
tion and wave amplitude, respectively. As might 
be expected, the greatest Doppler shifts occur at 
times closely adjacent to (but not precisely at) the 
moment the reflecting facet moves most rapidly 
across the path. It is interesting to note, however, 
that the rate of change of elevation and the rate of 
change of azimuth, corresponding to the first and 
second terms in the brackets of (8), oppose each 
other during this part of the cycle. For comparison, 
the two competing terms are included in figure 4b. 
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Figube 4. Angle of arrival and Doppler shift as functions of time for wavy layers with a single reflecting point. 



Presumably, slight departures from the idealized 
model assumed here could permit either term to 
predominate. Hence, rapid changes in Doppler at 
these moments are possible. 

5. Several Layers With a Single 
Reflecting Point 

When more than one layer is present, the curves 
for received signal strength are no longer as well- 
behaved as in figure 3. In figures 5 and 6 we have 
assumed two wavy layers — and in figures 7 and 8, 
three wavy layers — each supporting a single ray. In 
all of these cases we have assumed that the wavy 
layers are identical in shape but are at different 
elevations. On each figure we have shown not only 
the signal strength, but also the amplitude distri- 
bution and the power spectrum. 

Plots of typical signal strength (electric field) ver- 
sus time are shown in figures 5a to 8a. The rapid 
variations observed in these figures result from the 
beating between the two or three received rays. The 
beating results from the change in path lengths asso- 
ciated with different azimuth positions of the reflected 
ray. In figure 5a the two layers are at nearly the 
same height (around 1 deg above the horizon), and 
consequently both reflected rays are nearly equal in 
amplitude so that the interference fades are deep and 
relatively slow. However, in figure 6a the layers 
are farther apart (one at deg, the other at 0.75 
deg elevation) so that the reflection closer to grazing 
is appreciably stronger; the interference fades are 
therefore less pronounced, but the fading rate is 



markedly increased. In figures 7a and 8a, the 
addition of a third layer has made the received 
signal more erratic. The only difference between 
these two figures is the position of the middle layer. 
The fades are deeper if the middle layer is close to 
the lower layer (fig. 8), since the lower one has the 
strongest amplitude; but when the middle layer is 
close to the top one, the fading rate is seen to be 
more rapid. 

In figures 5b to 8b the signal distribution is plotted 
on a Rayleigh scale. On this scale a Rayleigh dis- 
tribution has a negative 45-deg slope as shown by 
the dashed curve. In figures 5 and 8 the distribu- 
tions roughly approximate a Rayleigh distribution — 
even though the signal is made up of two and three 
components respectively — and the phase relations 
are deterministic rather than random. In figures 7 
the distribution is a close approximation to a Rice 
distribution [Rice, 1945; Beckmann, 1961] which is 
the sum of a constant vector plus a Rayleigh distrib- 
uted vector. As has been noted previously for other 
circumstances, one finds that the distribution of a 
signal resulting from a small number of components 
can often approximate more complicated distribu- 
tions. 

In figures 5c to 8c we have plotted the power 
spectrum. The spectrum was obtained by com- 
puting the time autocorrelation and then taking the 
Fourier transform to find the power spectrum. A 
hanning window was used [Blackmail and Tukey, 
1958]. Though in most cases the power spectrum 
peaks up at only one or two frequencies, in figure 7 
quite a broad spectrum is observed. 
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Figure 5. Properties of received signal reflected 
from two layers that are close together. 
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Figure 6. Properties of received signal reflected 
from two layers that are far apart. 
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Figure 7. Properties of received signal reflected 
from three layers (two with high elevation). 



Figure 8, Properties of received signal reflected 
from three layers (two with low elevation). 
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6. Wavy Layers With More Than One 
Reflecting Point 

In figure 2 there are two cases where a single layer 
supports more than one reflected ray during at least 
part of the time (ry>l). Just as in the case of multi- 
ple layers, there is a beat between the several rays. 
An example of this effect is seen in figure 9 for the 
layer shown in the second case of figure 2. Only a 
detailed portion of the time record is shown, starting 
just before the wave moves from the position in 
which it supports one reflection to the position in 
which it supports three. With the onset of three 
reflections, rapid fading commences. It is faster — 
because of the rapid motions of the reflecting facets 
during this part of the cycle — than the fading associ- 
ated with multiple, relatively flat 0?<O) layers. 
However, this rapid fading occurs only when there 
are several rays present. The rest of the time the 
signal is quite steady. As mentioned earlier, the 
maximum slopes of the wavy layers must exceed a 
critical value before such rapid fading will occur. 




Figure 9. Rapid fading due to three rays received from a 
single layer. 

7. Conclusion 

We have examined quantitatively the conse- 
quences of an atmospheric model consisting of low- 
amplitude waves on a nearly horizontal surface of 
discontinuity — or abrupt change — in refractive index. 
The model is assumed bo be one deserving consid- 
eration as a mechanism for transhorizon propagation, 
and the results derived relate to phenomena which 
may be observed in transhorizon experiments. 
Specifically, we have evaluated the azimuth of 
arrival and its time rate of change, the signal level 
and its variation, the Doppler shift, and the received - 
signal distribution and power spectrum. We have 
investigated so far a variety of single-layer and 
multiple-layer conditions. 

The azimuthal deviations from a great-circle 
trajectory are strongly dependent on layer elevation 
as well as wave slope on the layer. Rate of change 
of azimuth is critically dependent on a parameter 



rj which is proportional to maximum wave slope and 
to layer elevation. For layers such that rj > 0.5, 
the velocity of the reflecting point on the layer may 
exceed the wave velocity; and for t?>1, multiple 
reflections from one layer may exist and their motions 
may have momentarily infinite velocities. Signal 
levels are influenced not only by the magnitude of 
the discontinuity and the grazing angle (in a manner 
very similar to that predicted by turbulent scat- 
tering), but also by layer curvature imparted by the 
wave. The consequent focusing results in signal 
peaks (of several decibels) coincident with the 
moments of fastest azimuth variation. Doppler 
shifts of a few cycles per second are likely — at 3 Gc/s 
for a 100-mile path — and are proportional to radio- 
frequency and to the first or second power of path 
length. These shifts change most rapidly at the 
moments of fast azimuth motion and signal peaking. 
The presence of more than one layer may lead to 
rapid fading: when two or more of the strongest 
reflections are nearly equal, the resulting signal-level 
distribution is approximately Rayleigh; and if one 
reflection predominates, it resembles a Rice distri- 
bution. Under some circumstances, particularly for 
a layer whose maximum slope just exceeds the 
critical value (?7>1), the fading may change abruptly, 
jumping to much higher rates during a small portion 
of the cycle. It is interesting to note that this type 
of phenomenon is similar to that attributed to the 
passage of an airplane across the transmission path. 
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